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Abstract 

Background: Huanglongbing (HLB) or citrus greening is a devastating disease of citrus. The gram-negative bacterium 
Candidatus Liberibacter asiaticus (Las) belonging to the a-proteobacteria is responsible for HLB in North America as well 
as in Asia. Currently, there is no cure for this disease. Early detection and quarantine of Las-infected trees are important 
management strategies used to prevent HLB from invading HLB-free citrus producing regions. Quantitative real-time 
PCR (qRT-PCR) based molecular diagnostic assays have been routinely used in the detection and diagnosis of Las. The 
oligonucleotide primer pairs based on conserved genes or regions, which include 16S rDNA and the (3-operon, have 
been widely employed in the detection of Las by qRT-PCR. The availability of whole genome sequence of Las 
now allows the design of primers beyond the conserved regions for the detection of Las explicitly. 

Results: We took a complimentary approach by systematically screening the genes in a genome-wide fashion, 
to identify the unique signatures that are only present in Las by an exhaustive sequence based similarity search 
against the nucleotide sequence database. Our search resulted in 34 probable unique signatures. Furthermore, 
by designing the primer pair specific to the identified signatures, we showed that most of our primer sets are 
able to detect Las from the infected plant and psyllid materials collected from the USA and China by qRT-PCR. 
Overall, 18 primer pairs of the 34 are found to be highly specific to Las with no cross reactivity to the closely related 
species Co. L. americanus (Lam) and Ca. L. africanus (Laf). 

Conclusions: We have designed qRT-PCR primers based on Las specific genes. Among them, 18 are suitable for the 
detection of Las from Las-infected plant and psyllid samples. The repertoire of primers that we have developed and 
characterized in this study enhanced the qRT-PCR based molecular diagnosis of HLB. 

Keywords: Detection system, Diagnostic, Candidatus Liberibacter asiaticus, Greening, Huanglongbing, Bacteria, Psyllid, 
Citrus 



Background 

Huanglongbing (HLB) or citrus greening is the most 
devastating disease of citrus, threatening the citrus in- 
dustry worldwide, and leading to massive reduction in 
fruit production as well as death of infected trees [1]. 
The causal agents of HLB are three closely related gram- 
negative, phloem-limited a-proteobacteria Candidatus 
Liberibacter species [2,3]. The heat tolerant strain Ca, L. 
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asiaticus (Las) is the most widespread in Asia as well as 
in the USA whereas Ca. L. americanus (Lam) is mostly 
limited to South America [2-4], Ca. L. africanus (Laf) is 
heat sensitive and localized to the African continent. All 
the three Liberibacter species are currently uncultured 
and are known to reside in the sieve tubes of the plant 
phloem [5] or in the gut of the phloem-feeding psyl- 
lids [6]. Psyllids are the natural vectors in transmit- 
ting the bacteria between plants [1,6]. The Asian 
psyllid, Diaphorina citri Kuwayama (Homoptera: Psyllidae) 
is responsible for transmitting Las and Lam in Asia and 
America, while the African citrus psyllid, Trioza erytreae 
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Del Guercio (Homoptera: Psyllidae), is the natural vector 
of Laf in Africa [7]. The characteristic symptoms of the in- 
fected plants include the yellow shoots, foliar blotchy mot- 
tles, along with poor flowering and stunting [1]. HLB also 
results in poorly colored, unpleasant tasting, reduced size 
fruit that shows staining of vascular columella and seed 
abortion [1]. Generally the fruit may remain partially green, 
for this reason HLB is also called citrus greening [1]. 
Chronically infected trees are sparsely foliated and display 
extensive twig or limb die-back and eventually die within 
three to five years [1]. Moreover, the disorders induced in 
diseased plants vary with cultivar, tree maturity, time of in- 
fection, stages of disease and other abiotic or biotic agents 
that affect the tree [1]. HLB symptoms also share certain 
similarities to nutrient deficiency [1], citrus stubborn dis- 
ease caused by Spiroplasma citri [8] and a HLB-like disease 
caused by a phytoplasma [9,10]. Early diagnosis and differ- 
entiation of Las infections from those defects and agents 
mentioned above, is thus critical to reducing the spread 
and devastation of this disease locally and via international 
trade, as well as minimizing the economic impact of poten- 
tial false positive diagnoses. 

Importantly, HLB and the Asian citrus psyllid (D. citri) 
are expanding to new citrus production areas. Currently, 
Asian citrus psyllid has been found in Florida, Texas, 
California, Arizona, Hawaii, Louisiana, Georgia, and 
Alabama in the USA, as well as in parts of South and 
Central America, Mexico, and the Caribbean. Meanwhile, 
HLB has not only been identified in Florida, Louisiana, 
South Carolina, Louisiana, Georgia, Texas and California 
of the USA; it has also been discovered in Cuba, Belize, 
Jamaica, Mexico, and other countries in the Caribbean 
[11]. While HLB and D. citri have been found in different 
producing areas, the number of infected trees and the 
psyllid vector population vary dramatically among dif- 
ferent regions. Thus, different strategies of manage- 
ment of HLB are recommended for different regions, 
according to the corresponding severity of HLB and 
occurrence of psyllid vectors. 

Currently, no efficient management strategy is avail- 
able to control HLB. For the recently Las-infected citrus 
producing areas such as California, prevention and 
eradication of HLB are the most efficient and cost- 
effective approaches. Additionally, Las infected trees are 
most often found to be asymptomatic during the early 
stage of infection. Thus, accurate early detection of Las 
in citrus plants and psyllids is critical for enacting con- 
tainment measures in non-endemic citrus producing 
areas. For the citrus producing areas without HLB, 
such as the Mediterranean region, accurate detection 
is critical for the success of quarantine measures 
against Ca. Liberibacter. 

Methods such as biological indexing using graft, dodder 
transmission [12], isothermal loop amplification (LAMP) 



[13], electron microscopy [1], DNA probes [14], enzyme- 
linked immunosorbent assays (ELISA) [15], conventional 
PCR [16-22] and quantitative real-time PCR (qRT-PCR) 
[22-26] are used for the diagnosis and confirmation of 
HLB. Although diagnostic tools like conventional PCR 
and LAMP showed good sensitivity, they were not consist- 
ent in detection of Las bacterium from infected plant and 
psyllid materials [6,13,25]. The current HLB diagnostic de- 
tection mainly employs qRT-PCR based methods due to 
their sensitive and quantitative nature. The initial qRT- 
PCR oligonucleotide primer sets for the detection of Las, 
targeted rplKAJL-rpoBC operon (p-operon: CQULA04f/r) 
[26], 16S ribosomal RNA gene (rDNA) (HLBasf/r) [23], 
EUB338f/EUB518r [27], ALF518f/ EUB518r [27] or spe- 
cies specific variable regions. EUB338f/EUB518r primers 
are universal to Eubacteria [27], while ALF518f/EUB518r 
primers identify a-proteobacteria universally [27] includ- 
ing Las, therefore not specific. Furthermore, the primers 
based on the conserved 16S and p-operon regions are 
popular but nevertheless have been shown to pose a po- 
tential specificity issue, as both false negatives and false 
positives have been reported [28]. Therefore, efforts have 
been directed towards developing effective qRT-PCR 
primers that target other non-conserved sequences. Re- 
cent studies made use of intragenic repeat regions of the 
prophage sequence for the detection of Las by qRT-PCR 
[25]. However, the intragenic repeat regions of the pro- 
phage sequence were also identified in Lam. Therefore, 
these primer pairs, hyvi/hyvii did not distinguish between 
Las and Lam, posing a specificity issue [25]. Consequently, 
primer pairs that specifically detect Las and make clear 
distinction among other phylogenetically closely related 
bacteria are essential. 

Here we took a complimentary approach to identify 
the genes that are unique to Las by a bioinformatic ana- 
lysis with the goal of expanding the arsenal of tools for 
Las detection. The advancement in the genome sequen- 
cing of Las [29] provides an opportunity to design 
primers based on species specific sequences for the de- 
tection of Las. We designed the oligonucleotide primer 
pairs specific to the identified unique genie signatures. 
We further validated their specificities and selectivity 
against closely related strains that demonstrated the ap- 
plication to Las-infected tissues and insect vectors by a 
qRT-PCR. 

Results and discussion 

Recently, the whole genome sequences of Las [29,30] 
have been sequenced. This allows for systematic screen- 
ing of unique Las genes in a genome-wide fashion. 
The availability of the genome sequences of the closely 
related species Lam [31], L. crescens (Lcr) [32] and Ca. 
L. solanacearum (Lso) [33], further effectively helps in 
identification of unique regions, by minimizing the 
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cross-species reactions, thereby enhancing the diag- 
nostic identification of Las in a more distinct manner. 

Bioinformatic analysis 

Several high-throughput applications have been devel- 
oped recently to design diagnostic primers using the 
whole genome sequence information including KPATH, 
Insignia, TOFI, and TOPSI [34-40]. Among them, KPATH, 
Insignia, and TOPSI have the potential to be used for de- 
sign of real-time PCR primers for qRT-PCR based assays 
for Las, whereas TOFI is used to design signatures for 
microarray-based assays. These methods mentioned above 
can be basically categorized into alignment-free and 
alignment-based approaches. The alignment-free approach 
uses both coding and non-coding regions of the genome 
and is useful for the genomes with less accurate sequence 
information, but generally result in high false positive rates 
as it does not involve pre-screening of the selected gen- 
omic loci for their discriminatory ability [37]. The 
alignment-based approach involves pre-screening of 
the selected genomic loci for their discriminatory ability 
[34]. This approach does not consider the genome annota- 
tion of genie and non-genic information, but rather aligns 
bigger regions of the genome, hence prone to lose shorter 
discriminatory sequence regions. Additionally, discrimin- 
atory ability of the selected regions are screened bioinfor- 
matically only on limited number of closely related 
species, which provide more opportunities for false posi- 
tives. We therefore took a complementary bioinformatics 
approach by pre-screening shorter genie regions against 
the nucleotide sequence database (nt) at NCBI, to identify 
all the possible unique genie regions from the Las gen- 
ome. The natural selection acts more strongly on genie re- 
gion, hence use of discriminatory sequences in this region 
results in less false positives as the organisms are under 
selection pressure [41]. Additionally, pre-screening against 
the nt is more effective as it contains the largest pool of 
well-annotated nucleotide sequences from different organ- 
isms. We envisioned that these two steps would result in 
more specific detection of target organism with less 
false positives, hence are included in our bioinformat- 
ics approach. 

There are -1100 genes assigned to the Las genome. 
Therefore, manual searching of each of these sequences 
against the nt database using BLAST program [42,43] is 
a laborious and time consuming procedure. Hence, we 
automated this sequence similarity search step by devel- 
oping a standalone PERL script (Additional file 1). This 
script performed the similarity searches for each of the 
Las gene against the specified database with hard-coded 
parameters for the BLAST program. Further, manual 
analysis of the resulting BLAST search output files is 
also laborious and time consuming; we therefore, auto- 
mated this step by developing a second PERL script 



(Additional file 2). This script automatically parsed all 
the BLAST output files and returned the Las sequences 
for which, no hits were found in other organisms. We 
refer to these sequences as probable unique sequences, 
because there are nearly no identical sequences found in 
other organisms (Figure 1). 

We performed the sequence similarity searches first by 
using stringent E- value of < 1 x 10~ 3 against nt database 
(Figure 1). This search resulted in -200 sequences that 
are unique to Las. This set of sequences is relatively high 
to validate experimentally; therefore, to further reduce 
the number of unique sequences, we performed the 
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Figure 1 Pictorial representation of the bioinformatics strategy 
employed to churn out the unique genie regions from Las 
genome. The input and output of each step are shown in oval or 
square boxes. Actions taken are noted to the left side of the arrow 
mark, while the information used is indicated to the right side of 
the arrow. 
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second sequence similarity search with a relaxed E-value 
of < 1. This search resulted in 38 unique sequences. The 
E-value of < 1 excludes the sequences with even little 
similarity to other organisms. Therefore, the resulting 38 
unique sequences are considered unique to Las and con- 
stitute the promising candidates for qRT-PCR based de- 
tection (Figure 1). 

We further searched the 38 unique sequences of Las 
against the phylogenetically closely related Lso, Lam, 
and Lcr. Because these organisms are closely related, we 
used the stringent E-value threshold of < 1 x 10" 3 for this 
similarity search. In order to achieve this E-value, the se- 
quences need to be highly similar between the Las, Lso, 
Lam, and Lcr. Therefore, this close species filter proced- 
ure potentially eliminates all the Las sequence targets 
that could lead to false positive results in qRT-PCR 
based molecular diagnostic assays. Consequently, we fur- 
ther eliminated four conserved sequences from the list 
of 38 unique sequences, resulting in a total of 34 poten- 
tial sequence signatures. We could not apply this close 
species filter step against Laf genome as its genome is 
yet to be sequenced. 

Five (-15%) of the 34 unique gene sequences namely 
CLIBASIA_05545, CLIBASIA_05555, CLIBASIA_05560, 
CLIBASIA_05575 and CLIBASIA_05605 are in the pro- 
phage region of the Las genome. All these five unique 
sequences are located upstream of the genomic locus 
CLIBASIA_05610 encoding a phage terminase. There 
are possibly 30 genes that represent the complete pro- 
phage genome within the Las genome [25,44], of which 
16 open reading frames (ORFs) are upstream of the 
phage terminase, while the remaining 13 ORFs are down- 
stream. The prophage genes CLIBASIA_05610 (primer 
pair 766 F and 766R) and CLIBASIA_05538 (primer 
pair LJ900F and LJ900R) have been targeted in previous 
studies by both conventional as well as qRT-PCR based 
assays [25,44]. 

We further analyzed the genomic orientation of the 34 
unique genes. This analysis revealed that 15 (-44%) of 
them are oriented on the sense strand, while the remaining 
19 (-56%) were present on the anti-sense strand (Add- 
itional file 3: Figure SI). The sequence length of these 
unique genes ranged from 93 to 2595 base pairs (bp) 
(Additional file 4: Table SI). 

Designing of Las specific primers and experimental 
validation of the specificity and sensitivity of qRT-PCR 
assay to detect Las 

Based on the genome sequence of Las strain psy62, we 
designed 34 qRT-PCR primer pairs that specifically tar- 
get the 34 unique sequences identified in our bioinfor- 
matic analyses (Additional file 4: Table SI). We designed 
the melting temperature (Tm) of each of these primers 
to range from 59°C to 65°C with an optimum of 62°C. 



The GC content of the primers ranged from 35% to 65% 
with an optimum of 50%. The PCR amplicon sizes for 
each primer set are between 84 to 185 bp (Additional 
file 4: Table SI). 

In addition to the novel primers designed in this work, 
we also used a set of control primers that have been pre- 
viously used in a qRT-PCR based detection of Las. These 
known primers include 16S rDNA pairs specific to the 
three different Candidatus Liberibacter species (HLBasf/ 
r: Las, HLBamf/r: Lam and HLBaf/r: Laf) [23], p-operon 
(CQULA04f/r: p-operon) [26], intragenic repeats regions 
of the prophage sequence (LJ900f/r: Prophage) [25], and 
the primer pair specific to the plant cytochrome oxidase 
(COXf/r: COX) gene [23] as a positive endogenous 
control. 

We performed qRT-PCR assays to test the specificity 
of the designed primers using total DNA extracted from 
Las-infected citrus plants as a template. To further val- 
idate the specificity of these primers, we also included 
total DNA from the phylogenetically closely related spe- 
cies Lam and Laf in our test. Additionally, DNA ex- 
tracted from healthy citrus plant was used as a negative 
control, whereas water served as a no template control. 
The results of qRT-PCR assays are listed in Table 1. 

Most of our novel custom designed primer pairs tar- 
geting the unique gene sequences were indeed found to 
be highly specific to Las, as assessed by qRT-PCR assays 
(Table 1). Among the 34 primer pairs, 29 produced 
amplicons only when Las-infected citrus plant DNA was 
used as a template, with an average C T value ranged 
from 19.48 to 27.47. Two primer pairs, P13 and P15, 
didn't produce any amplicons under the standard condi- 
tions tested. The other three primer pairs, PI 9, P27 and 
P28, produced amplicons when Las or Laf infected plant 
DNA was used as a template, indicating PI 9, P27 and 
P28 could be used to detect both Las and Laf. We were 
unable to filter for cross-reactivity of PI 9, P27 and P28 
in the bioinformatic analysis, because the Laf genome 
sequence is currently unavailable. With the exception of 
these three primer sets that showed amplicons with Laf 
template, none of the other primer sets produced any 
amplicons with DNA of Lam, Laf, and healthy citrus or 
water as template, which further confirms the specificity 
of these primers to the Las. 

We further evaluated the specificity of these primer 
sets using DNA templates from various citrus associated 
fungal and bacterial pathogens including Colletotrichum 
acutatum KLA-207, Elsinoe fawcettii, Xanthomonas axo- 
nopodis pv. citrumelo 1381, X. citri subsp. citri strains 
306, A w , and A*. Only two primers sets, P20 and P21 
showed unspecific amplification against template DNA 
extracted from fungal pathogen C. acutatum KLA-207 
(Table 1). C. acutatum causes citrus blossom blight, 
post-bloom fruit drop and anthracnose symptoms that 
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Table 1 Specificity and sensitivity of the novel primers in the detection of Las as shown by qRT-PCR assay 

Primer Target gene Las C T value of the qRT-PCR # 



Negative control Other controls 







C T value 


R 2 value + 


Slope + 


Laf 


Lam 


Healthy plant 
tissue 


Water 


CI 


C2 


C3 


C4 


C5 


C6 


PI 


CLIBASIA_05555 


20.54 


0.9944 


-0.2883 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P2 


CLIBASIA_04315 


19.99 


0.9867 


-0.2849 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P3 


CLIBASIA_05575 


20.15 


0.9991 


-0.2847 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P4 


CLIBASIA_05465 


19.52 


0.9618 


-0.2897 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P5 


CLIBASIA_01460 


19.48 


0.9995 


-0.2969 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P6 


CLIBASIA_05145 


22.29 


0.9971 


-0.3057 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P7 


CLIBASIA_05545 


20.11 


0.9972 


-0.3407 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P8 


CLIBASIA_05560 


19.92 


0.9982 


-0.3132 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P9 


CLIBASIA_02025 


20.12 


0.9875 


-0.2743 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P10 


CLIBASIA_05605 


20.18 


0.9945 


-0.2781 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P11 


CLIBASIA_03090 


23.61 


0.9997 


-0.2867 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P12 


CLIBASIA_03875 


27.47 


0.9992 


-0.2563 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P13 


CLIBASIA_02305 


UD 


NT 


NT 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P14 


CLIBASIA_05495 


21.25 


0.9974 


-0.2594 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P15 


CLIBASIA_02660 


UD 


NT 


NT 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P16 


CLIBASIA_02715 


20.26 


0.941 1 


-0.3480 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P17 


CLIBASIA_031 10 


20.11 


0.9994 


-0.2786 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P18 


CLIBASIA_03675 


20.02 


0.9967 


-0.2780 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P19 


CLIBASIA_03725 


19.91 


NT 


NT 


35.29 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P20 


CLIBASIA_03955 


21.08 


NT 


NT 


UD 


UD 


UD 


UD 


37.41 


UD 


UD 


UD 


UD 


UD 


P21 


CLIBASIA_04030 


20.30 


NT 


NT 


UD 


UD 


UD 


UD 


32.93 


UD 


UD 


UD 


UD 


UD 


P22 


CLIBASIA_04150 


24.00 


NT 


NT 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P23 


CLIBASIA_04310 


20.76 


0.991 


-0.2976 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P24 


CLIBASIA_04330 


20.85 


0.9986 


-0.2635 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P25 


CLIBASIA_04405 


21.60 


0.9987 


-0.3051 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P26 


CLIBASIA_04425 


20.41 


0.9994 


-0.3032 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P27 


CLIBASIA_02645 


21.77 


NT 


NT 


38.61 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P28 


CLIBASIA_04515 


22.00 


NT 


NT 


38.63 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P29 


CLIBASIA_04530 


19.00 


0.9919 


-0.2852 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P30 


CLIBASIA_04550 


22.48 


0.9938 


-0.2708 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P31 


CLIBASIA_05230 


21.68 


0.9941 


-0.2771 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P32 


CLIBASIA_05480 


21.48 


0.988 


-0.2776 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


UD 


P33 


CLIBASIA_04475 


20.84 


0.9913 


-0.2644 


UD 


UD 


UD 


UD 


NT 


UD 


UD 


UD 


UD 


UD 


P34 


CLIBASIA_05505 


22.70 


0.9893 


-0.2791 


UD 


UD 


UD 


UD 


NT 


UD 


UD 


UD 


UD 


UD 


CQULA04F/R 


(3-operon 


22.11 


NT 


NT 


UD 


UD 


UD 


UD 


NT 


NT 


NT 


NT 


NT 


NT 


LJ900f/r 


Prophage 


22.25 


NT 


NT 


UD 


UD 


UD 


UD 


NT 


NT 


NT 


NT 


NT 


NT 


HLBas/r 


16Sas 


24.33 


0.9998 


-0.3057 


NT 


NT 


UD 


UD 


NT 


NT 


NT 


NT 


NT 


NT 


HLBam/r 


16Sam 


NT 


NT 


NT 


NT 


24.68 


UD 


UD 


NT 


NT 


NT 


NT 


NT 


NT 


HLBaf/r 


16Saf 


NT 


NT 


NT 


21.28 


NT 


UD 


UD 


NT 


NT 


NT 


NT 


NT 


NT 


COXf/r 


Cox 


14.80 


NT 


NT 


15.21 


18.54 


16.15 


UD 


NT 


NT 


NT 


NT 


NT 


NT 



+ Las-infected psyllids total DNA was serially diluted spanning up to five logs and used as a template in the qRT-PCR assay. R 2 and slope were further calculated 
from a plot of C T values versus log dilution factor. # qRT-PCR was conducted by using template DNA samples of Las, Laf, Lam, CI: Colletotrichum acutatum KLA-207, 
C2: Elsinoe fawcettii, C3: Xanthomonas axonopodis pv. c/frame/01381 , C4: Xanthomonas citri subsp. citri A w , C5: Xanthomonas citri subsp. citri A*, C6: Xanthomonas 
citri subsp. citri 306. The C T values are average of three replicates for each primer pair. UD: undetected; NT: Not tested. 



Table 2 qRT-PCR detection of Las from plant samples that were collected from different locations in USA and China 



Primer 
pairs 












C T value of qRT-PCR using infected plant DNA samples as template* 




















DNA samples from Florida, USA 












DNA samples from China 




Home stead 


Orange 


Polk 


Lake wales 


1 1 • LI- 1 — 

Highlands 


de Soto 


St Lucie 


Hendry 


Hickory 


Hardee 


Charlotte 


Indian river 


Hai nan 


Jiang xi 


Guang xi 
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Guang dong 


P1 


23.46 


22.24 


25.33 
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z4./ z 
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ZD.OJ 
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oa nn 

ZD.UU 
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i in 
UU 


QO 

32.DD 
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UD 


P2 
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23.07 


oa ao 

zo.yu 


OQ Q 1 


iq qo 
zj.oU 


OQ 07 
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UU 
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Q 1 /I Q 
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07 Q1 
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P5 


24.44 
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22.47 
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ZD.U/ 
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ZO. I / 
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z4.4j 


OQ AO 
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28.91 


28.53 


26.66 
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27.31 


O^ OO 
ZJ.UZ 


Q1 AQ 
O I .Do 


OQ AQ 

zo.4y 


26.98 


P6 
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28.02 
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zy.uo 


oc; 07 
2d. 2/ 
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P7 
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UD 
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ZD.UU 
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UU 


Q 1 1 Q 

o i . i y 
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Z/.4U 


UD 
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25.96 
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UU 


Q 1 7/I 
O I ./4 


QO Q7 

ou.y/ 


UD 
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25.89 


24.02 


28.51 


24.84 


OQ 

ZO.JJ 


QO ^O 
OU.JZ 


OA AH 


QQ t^O 

oU.jz 


31 .72 


30.66 


28.08 


30.54 


28.47 


OA OQ 

zD.uy 


Q7 £^A 
3/ .DO 


Q ^ A 1 
DD A I 


29.28 


P16 


25.50 


23.36 


27.87 


23.20 


OA QC 
ZO.OJ 


OQ /I 1 

zo.4 I 


ZD.D/ 


OO 1 Q 

zy. I O 


29.41 


29.54 


27.57 


28.88 


28.10 


nr QO 
ZJ.OZ 


QO QA 


07 07 

21 .21 


27.81 


P17 


25.95 


24.09 


28.18 


23.65 


Z/.D4 


OQ QA 


OA A1 
ZD.D I 


OQ QQ 

zy.yu 


29.50 


31.09 


28.14 


30.92 


29.34 


07 Q1 

z/.U I 


QA 1 O 
oD. I Z 


QO OQ 

oU.zo 


29.20 


P18 


25.17 


23.1 1 


28.02 


23.07 


2/ A3 


OQ 7^ 
2(3./ D 


Oc; QQ 

zj.yy 


OQ QA 

zo.yD 


29.36 


29.15 


28.19 


29.09 


28.67 


OA A 1 

zD.4 I 


QO 1 7 
OZ. I / 


07 QQ 

z/.oy 


28.79 


P23 


26.41 


24.05 


29.28 


24.35 


zo.U4 


qo oo 
oU.zz 


27 '.75 


31.15 


32.14 


32.95 


29.77 


31 .48 


30.31 


07 A7 
Z/.D/ 


QA 7Q 

oD./o 


QO QA 

oU.oD 


30.63 


P24 


26.14 


23.83 


28.80 


23.68 


27 58 


29 68 


27 28 


30 86 


32.14 


31.87 


30.71 


31.84 


29.75 


27 51 


37 70 


30 80 


30.05 


P25 


25.04 


22.68 


27.97 


22.90 


26.67 


28.28 


25.92 


28.63 


OQ O/l 

zy.U4 


QO QO 

oU.oU 


07 QO 

z/.oU 


OQ 77 

zy./ / 


27.81 


25.47 


36.49 


29.31 


29.31 


P26 


25.1 1 


23.1 1 


27.65 


22.86 


27.31 


28.53 


25.71 


28.55 


29.57 


28.66 


27.89 


29.49 


28.41 


26.20 


31.67 


27.50 


28.38 


P29 


24.73 


22.72 


27.21 


22.60 


26.65 


27.85 


25.42 


29.36 


29.56 


29.28 


27.17 


29.13 


27.39 


25.33 


34.12 


28.03 


27.51 


P30 


26.46 


24.87 


30.59 


24.55 


28.91 


30.73 


27.79 


29.69 


31.25 


31.89 


28.33 


30.69 


29.32 


26.60 


35.91 


29.90 


30.71 


P31 


27.19 


25.05 


29.83 


24.77 


29.43 


31.03 


27.88 


31.23 


32.67 


31.14 


29.94 


30.71 


30.28 


27.96 


34.28 


29.94 


31.58 


P32 


26.65 


24.65 


29.13 


23.73 


28.24 


29.40 


25.93 


29.44 


30.58 


30.20 


28.11 


29.82 


28.94 


26.60 


33.83 


29.23 


28.77 


P33 


25.55 


23.35 


28.08 


23.33 


27.03 


28.42 


26.32 


30.32 


30.58 


30.36 


27.83 


29.79 


28.41 


25.80 


32.99 


30.71 


28.37 


P34 


26.49 


24.29 


29.62 


24.46 


28.14 


29.45 


26.22 


28.50 


29.66 


30.85 


26.67 


29.28 


27.24 


25.66 


36.14 


29.07 


29.52 


HLBas/r 


24.76 


22.97 


27.55 


22.80 


31.02 


29.94 


27.24 


27.45 


28.02 


27.20 


28.90 


27.95 


27.06 


25.04 


30.40 


25.93 


25.78 



# Las-infected plant DNA samples were collected from 12 different locations in Florida, USA, and 5 different locations in China. The C T values indicated are average of three replicates for each primer pairs. 



UD: Under determined. 
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are phenotypically distinguishable from citrus HLB. The 
P20 and P21 were not filtered by the bioinformatic ana- 
lysis since C. acutatum genome sequence was unavail- 
able in the database. Because of the complexity of the 
natural microbial community and the limited number of 
sequences available in the current nucleotide sequence 
database, it is impossible to completely filter out all the 
potential false positives bioinformatically. However, false 
positives could be identified experimentally by combin- 
ing the different sets of primer pairs by a consensus ap- 
proach [37]. We eliminated these two primer sets from 
further evaluation in this study. 

The melting temperature analysis of the amplicons 
produced from our novel primer set with Las as a tem- 
plate indicated that amplicons were of a single species. 
This suggests that there is no off target amplification for 
our primer pairs on the Las genome. Overall, the experi- 
mental validation of the 34 novel primer sets specific to 
unique targets revealed that 27 (-80%) of these targets 
are indeed specific to the Las genome (Table 1). This 



demonstrates the significance of the bioinformatics strat- 
egy employed here for identifying the suitable target re- 
gions for the detection of the bacteria by qRT-PCR 
based methods. These 27 novel primer pairs were se- 
lected for further characterization. 

To test the sensitivity of our designed novel primers, 
serial dilutions of Las-infected psyllid DNA was used as 
a template in the qRT-PCR assay. This serial dilution 
qRT-PCR assay indicated that most of our novel primer 
pairs were able to detect Las up to 10 4 dilutions from 
the initial template DNA concentration, which is com- 
parable to that of the primer set targeting Las 16S 
rDNA (Table 1). However, lower sensitivity was ob- 
served in the case of primer pairs P9, P12, P14 and 
P22, which were eliminated from further study. The 
remaining 23 primer pairs were able to detect Las up 
to 10 4 dilutions, with a correlation co-efficient (R 2 >0.94) 
between the C T values and dilutions (Table 1). This dem- 
onstrates the high sensitivity of these 23 primers in the de- 
tection of Las. 




©Orange 

©Orange 

©Polk 

©Polk 

©CREC 

©Lake wales 

©Hickory 

©Indian River 

©Hardee 

®De Soto 

©Charlotte 

10a Highlands 

©Highlands 

@ St. Lucie 

©Hendry 

©Miami 

(14) Homestead 



< ^\^ c Las from psyllids 

Las from plant samples 



Not on scale 



Figure 2 Schematic representation of the plant and the psyllid samples collected from Florida. Las-infected plant DNA samples were 
collected from 12 different locations and psyllids from 5 different locations in Florida, USA. The color shaded symbols for representative plant and 
psyllid samples are based on their average infection level across all the primer pairs tested based on C T values. 
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qRT-PCR detection of Las from plant and psyllid DNA 
samples isolated from diverse locations in USA and 
China 

In order to further demonstrate the degree of applicabil- 
ity of the 23 primer pairs in the detection of Las from 
infected biological material, we performed qRT-PCR on 
the various Las-infected plant and psyllid DNA samples. 
Considering the potential variation in nucleotide se- 
quences of Las isolates in different geographic locations 
that might affect our detection due to the potential nu- 
cleotides changes of the target unique genes, we col- 
lected Las-infected plant DNA samples as tabulated in 
Table 2, from not only USA, but also from China, where 
Las was reported more than 100 years ago [1]. We tested 
the 23 primer pairs on 17 Las-infected plant DNA sam- 
ples. Of these 17, 12 were collected from different loca- 
tions in Florida, USA (Figure 2, Table 2), and the 
remaining five were collected from different locations in 
China (Table 2). Additionally, Las-infected psyllid DNA 
samples collected from five different locations in Florida, 
USA, were also included in the qRT-PCR assays (Table 3, 
Figure 2). 

All the 23 primer pairs detected Las from all 12 Flor- 
ida HLB diseased plant samples (Table 2) and 5 psyllid 
DNA samples (Table 3) in a qRT-PCR assay, which fur- 
ther validated the detection applicability of our novel 
primers (Figure 2). However, 4 of the 23 primer pairs 
(PI, P7, P8 and P10) failed to produce amplicons with the 
infected plant DNA sample from Jiangxi and Guangdong 
Province, China (Table 2). Primer pair P3 produced no 
amplicon with Jiangxi sample, and produced unspecific 
amplicon with the Guangdong sample (with an altered 
PCR product size, data not shown). Interestingly, all these 
5 primer pairs target the genes located in prophage region 
of the Las genome (Additional file 3). These primers (PI, 
P3, P7, P8 and P10) based on prophage genes could detect 
Las from Florida, but not from Jiangxi and Guangdong 
province, China. This is consistent with previous report 
[44], that prophage was detected in only 15.8% of the 120 
HLB diseased citrus samples acquired in Guangdong 
Province, China, but was detected in 97.4% of the 39 Las 
positive citrus samples acquired in Yunnan Province, 
China. This suggests that those prophage genes are not 
universally present in all strains of Las. Alternately, the 
prophage sequences were found to be highly variable 
among the strains tested. 

Conclusions 

We have successfully designed 18 novel primer pairs, 
which are specific to Las. These primers will provide an 
additional arsenal to qRT-PCR based detection of Las- 
infected plants and psyllids. Compared to the commonly 
used primers based on 16S rDNA and (3-operon, the 
18 primers developed in this study have comparable 



Table 3 qRT-PCR detection of Las from psyllid DNA 
samples that were collected from different locations in 
Florida, USA 



Primer C T value of qRT-PCR using infected psyllid DNA 

pairs samples as template* 
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P25 


33.15 


24.10 


28.46 


26.78 


24.77 


P26 


33.40 


25.59 


29.74 


28.07 


25.58 


P29 


33.42 


25.14 


29.49 


27.73 


25.29 


P30 


36.28 


26.53 


32.12 


29.65 


27.07 


P31 


36.10 


27.13 


31.67 


29.94 


27.43 


P32 


35.53 


26.40 


31.06 


29.22 


27.23 


P33 


33.86 


25.01 


30.00 


27.92 


25.65 


P34 


34.99 


25.74 


30.93 


28.58 


26.43 


HLBas/r 


33.41 


25.10 


29.09 


27.86 


25.57 



# Las-infected psyllid DNA samples were collected from 5 different locations in 
Florida, USA. The C T values indicated are average of three replicates for each 
primer pair. 



sensitivity. Moreover, these primers could successfully 
differentiate Las from Lam, Laf and other common 
microbes associated with citrus. 

Methods 

Bioinformatics 

The nucleotide sequences of Las with accession number 
NC 012985 [29,45], Lso with accession number NC 014774 
[33], Lcr with accession number NC_019907 and compre- 
hensive nucleotide (nt) database (26 th July 2012) were 
downloaded from the NCBI ftp server (ftp.ncbi.nih. 
gov). The stand-alone BLAST [42,43] was used to 
search the Las genes against nt, Lso and Lcr databases 
using a custom-made PERL script 1 (Additional file 1) 
by varying the E-value with all other parameters kept 
to a default value. The output files of the BLAST 
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searches were further parsed using a second custom- 
made PERL script 2 (Additional file 2). 

Plant and psyllid materials and extraction of DNA 

Las infected citrus leaf samples with typical visible 
symptoms were collected from 2 years old infected sweet 
orange (Citrus sinensis) plants maintained at the Citrus 
Research and Education Center (CREC), Lake Alfred, 
Florida, USA. As a negative control, the leaves from 
healthy citrus plants were collected from pathogen-free 
seedlings grown in the healthy plant greenhouse main- 
tained at CREC, Lake Alfred, Florida, USA. The Laf and 
Lam infected samples were obtained from South Africa 
and Brazil respectively. The total DNA from the leaves 
of citrus was extracted using the protocol mentioned 
elsewhere [46]. Briefly, the leaves were washed under tap 
water and surface sterilized in 35% bleach (2% active 
Chlorine) and 70% (v/v) ethanol for 2 min each. The 
sterilized leaves were further rinsed three times in sterile 
water. The midribs from the leaf samples were separated 
and cut into small pieces. Approximately 100 mg of 
midrib pieces were used from each sample to extract the 
DNA using the Wizard® genomics DNA purification kit 
(Promega, Madison, WI, USA). The extracted DNA was 
suspended in 100 ul H 2 0. 

Las infected psyllids {Diaphorina citri) were maintained 
on confirmed Las-infected sweet orange plants at the 
CREC, Lake Alfred, FL, USA. In this work, 16 psyllids 
(around 20 mg) were pooled and the total DNA was 
extracted using a DNeasy Blood & Tissue Kit (Qiagen, 
Valencia, CA). The extracted DNA was suspended in 
100 ul H 2 0. The quality and quantity of the extracted DNA 
was determined using a NanoDrop™ 1000 spectrophotom- 
eter (NanoDrop Technologies, Inc., Wilmington, DE). 

Quantitative real-time polymerase chain reaction (qRT-PCR) 

Gene specific primers were designed using Primer- 
QuestSM from Integrated DNA technologies (IDT), 
Coralville, Iowa (Additional file 4: Table SI). qRT-PCR 
experiments were performed using ABI PRISM 7500 
FAST Real-time PCR System (Applied Biosystems, Foster 
City, CA, US) in a 96-well plate by using an absolute 
quantification protocol. The reaction mixture in each well 
contained 12.5 uL 2x FAST SYBR® Green PCR Master 
Mix reagent (Applied Biosystems), 2 uL DNA template 
(-30 ng), 0.625 uL of 10 uM of each gene-specific primer 
pair in a final volume of 25 uL. The standard thermal pro- 
file for all amplifications was followed, which involved 95°C 
for 20 min followed by 40 cycles of 95 °C for 3 sec, and 
50°C for 30 sec. All assays were performed in triplicates. 

Melting curve analysis was performed using ABI PRISM 
7500 FAST Real-time PCR System Software version SDS 
vl.4 21 CFR Part 11 Module (Applied Biosystems") to 
characterize the amplicons produced in a PCR reaction. 



Additional files 



Additional file 1: PERL script 1 facilitates the similarity search in an 
automated fashion. This script performs similarity searches against the 
specified nucleotide sequence database using a stand-alone BLAST program 
for each of the input gene sequences from the Las genome. 

Additional file 2: PERL script 2 facilitates the identification of 
unique genes to Las. This script facilitates the identification of unique 
genes by automatically parsing all the BLAST output files generated from 
the Additional file 1 PERL script land returns the unique gene sequences 
with no similarity to the DNA sequences of other organisms. 

Additional file 3: Figure SI. Snapshot of the unique genes identified 
by bioinformatics is shown in the context of the whole genome of Las. 
The absolute positions of the regions are shown. The novel unique 
regions of Las identified in this study are shown in bluish green, while 
the currently known targets are colored in green. 

Additional file 4: Table SI. Custom designed primer pairs specific to the 
unique sequences of Las identified by bioinformatic analysis. The forward 
and reverse primer pair for each of the unique genie regions is given. The 
product size for each of the primers is shown along with the %GC content. 
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